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ABSTRACT 

We present a version of the Fourier Bessel method first introduced by [Fisher. 



Lahav et al (199¥)| and ^aroubi et al (1995)| with two extensions: (a) we amend 



the formahsm to allow a generic galaxy weight which can be constant rather 
than the more conventional overweighting of galaxies at high distances, and (b) 
we correct for the masked zones by extrapolation of Fourier Bessel modes rather 
than by cloning from the galaxy distribution in neighbouring regions. We test 
the procedure extensively on A^-body simulations and find that it gives generally 
unbiased results but that the reconstructed velocities tend to be overpredicted 
in high-density regions. Applying the formalism to the PSZz redshift catalog, 
we find that (3 = 0.7 ± 0.5 from a comparison of the reconstructed Local Group 
velocity to the CMB dipole. From an anisotropy test of the velocity field, 
we find that {3 = 1 CDM models models normalized to the current cluster 
abundance can be excluded with 90% confidence. The density and velocity 
fields reconstructed agree with the fields found by [Branchini et al (1998)| in most 
points. We find a back-infall into the Great Attractor region (Hydra- Cent aurus 
region) but tests suggest that this may be an artifact. We identify all the 
major clusters in our density field and confirm the existence of some previously 
identified possible ones. 



Subject headings: Cosmology: theory - galaxies: clustering, - large-scale 
structure, large-scale dynamics 
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1. Introduction 

Redshift surveys provide the only possibility for determining the three-dimensional 
density field of luminous matter which is crucial for studies of mass concentrations, the 
power spectrum, dynamical analyses to probe the relationship between dark and luminous 
matter and many other areas of observational cosmology. However, the relation between 
redshift space and real space, although given straightforwardly by Bubble's law in the limit 
of large distances, is distorted on smaller scales by galaxy peculiar velocities (a full analysis 
of this was first given by [Kaiser (1987)| ). Hence correcting for the distortions becomes the 
most important physical problem associated with redshift surveys. Several different types 
of these surveys are currently available, varying widely in their depth, sky coverage, and 
sampling density. Since this paper will be mainly concerned with analysis of the density and 
velocity fields, we use the recently completed PSC2; survey, which is ideal for our purposes 
because of its large number of galaxy redshift and its near-complete sky-coverage. Its depth 
is small enough to render effects of space curvature and galaxy evolution negligible in our 
calculations. 

There are several methods of correcting for distortions to recover the real space density 
and velocity fields. These all use either linear theory or the Zel'dovich approximation, and 
are therefore ultimately very limited in their ability to reconstruct the high-density regions.0 
They can roughly be separated into two types of methods: 

Iterative methods: Since peculiar velocities are caused by gravitational acceleration, 
the velocity field can be recovered from the density field. Iterative methods use the redshift 
space density field to calculate a peculiar velocity field which can then be used to correct 
the density field distortions. The procedure is repeated until the velocity field converges - 
see [Yahil, Strauss, et al (1991) and Kaiser, Efstathiou, Ellis et al (1991] for slightly different 



versions of this method. 

Basis function methods: When transforming the measured redshift space density field 
into a combination of angular and radial basis functions, the distortion is concentrated in 
the radial part and its correction becomes an algebraic matrix problem. There are several 
different versions of this approach: [Nusser and Davis (1991)| transform the angular part 



into basis functions but express the radial part in differential equations which they then 
solve numerically. Their method uses the Zel'dovich approximation. Pisher, Lahav et~ai 



(1994)1 (FLHLZ in the following) and [Zaroubi et al (1995)| transform both angular and 



radial parts into basis functions, using a combination of spherical harmonics and spherical 
Bessel functions. 



^However, this may change when techniques based on the fully nonlinear variational 



method of Peebles (1989) are developed to deal with large redshift surveys. The most recent 



work in this area is by [Sharpe et al (1999) . 
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Other approaches include [Kaiser and Stebbins (1991)| and [Saunders and Ballinger 



1998 



Work on nonlinear corrections for the evolution of the power spectrum includes [Peacock 



and Dodds (1994)] , [Fisher and Nusser (1996)] , and [Taylor and Hamilton (1996)[ ; a critical 



discussion of these may be found in [Hatton and Cole (1998 



In this paper, we extend and apply the Fourier Bessel method which was first 
introduced by FLHLZ and Zaroubi et al (1995) . Their treatment of dynamics using this 
expansion is new, but the idea of using spherical harmonics and even spherical Bessel 
functions goes back to Peebles (1973)[ It has been apphed to redshift surveys by ^charf and 
Scharf, Hofmann, Lahav, Lynden-Bell (1992)[ , and [Fisher, Scharf and Lahav 



Lahav (1993 



(1994)| . We will follow the original method closely but introduce two extensions: 



Conventionally, the weight given to each galaxy increases as the number density of 
galaxies decreases (as it does with radius in a flux-limited survey). We generalise the 
formalism to allow constant or any other weight, since we feel that care has to be 
taken with the conventional procedure: The purpose of this type of weighting is to 
exaggerate the mass of galaxies at higher radii where the sampling is poor, thereby 
making it possible to determine a density fleld. Note however, that since this also 
exaggerates the shot noise, the density fleld at high radii is subsequently considered 
unreliable and fluctuations are smoothed away. 

We correct for the mask by also using a basis function approach rather than the 
more usual cloning mechanism. The difference between these two methods is not 
an issue for IRAS-based surveys, which have very good sky coverage, but will be 
more interesting in the future with the advent of very deep surveys of small angular 
coverage such as the 2dF survey. 



For our analysis, we use the recent PSC2; redshift survey. It contains approximately 
15,500 galaxies (almost all) detected in the IRAS Point Source Catalog ( [Saunders (1996)| , 
Saunders et al (1998)[ ) with 60 /xm flux larger than 0.6 Jy. Our subsample contains 10549 
PSC2; objects within 170 /i^^Mpc and with positive galaxy identiflcation. Regions not 
surveyed by IRAS (two thin strips in ecliptic longitude and the area near the galactic plane 
deflned by a V-band extinction of > 1.5 mag) are excluded from the catalog which therefore 
covers ~ 84 % of the sky. 

The structure of the paper will be as follows: In section Q, we describe the correction 
for redshift space distortions using the Fourier-Bessel set of basis functions, and discuss 
the Wiener Filter smoothing procedure for suppressing shot noise. Details are given in the 
Appendices. In all of this, we mostly follow FLHLZ and [Zaroubi et al (1995)[ , but we have 
chosen to give derivations in full since we extended the original formalism and also changed 
some normalisations to render them more intuitive. In section |] we present our method of 
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correcting the masked areas. In section ^ we then test the method on A^-body simulations 
to evaluate errors and systematic biases. 

Sections ^, |^, and |^ present the analysis of the PSCz redshift survey. In section ^, we 
recover (3 from a comparison between the reconstructed velocity of the Local Group and 
its value known from the dipole anisotropy in the Cosmic Microwave Background (CMB). 
This procedure has a long tradition since the Local Group velocity is one of the few that 
is known accurately enough for a meaningful comparison to be made to its reconstructed 
value. However, the one-number-statistic nature of the procedure also makes it highly 
susceptible to systematic biases and we therefore take particular care to evaluate the error 
by analysing the scatter in the A^-body simulations. 

Section |] analyses the magnitude of the redshift space distortions to recover another 
estimate of /5 . It was first pointed out by [Kaiser (1987) that the distortions themselves 



obviously depend on /? so that by analysing their magnitude, it should be possible to 
recover a value for that parameter. Several versions of such analyses have been done to 
date - ours relies on the fact that, if we correct the redshift space distortions assuming a 
wrong value for (3 , the resulting density and velocity fields will be anisotropic, i.e., there 
will be a systematic difference between the radial and the other two directions. We develop 
a simple test for this anisotropy and estimate (3 for the PSC2; again by comparing to the 
corresponding results for A^-body simulations. 

In section |^ we discuss the reconstructed density and velocity fields respectively and 
compare them to other recent reconstructions using the PSC2; and similar catalogs. Both 
of these fields result naturally from the reconstruction method we use if in a somewhat 
smoothed form due to the fact that our formalism only works in the linear regime. 



2. The Fourier-Bessel method 

The idea of this method is to express the overdensity as a Fourier-Bessel expansion 

(5(r, w) = ^ yim{^)3i{hnr)5imn- (1) 

Imn 

Here us denotes the angular coordinates, Yim are spherical harmonics, ji are spherical Bessel 
functions, kin are a set of wavenumbers that depend on the boundary conditions assumed, 
and 5imn are the expansion coefficients. Once the expansion coefficients are known,|^. the 



^"Since the expansion has to be finite, the radial and angular resolutions are finite. The 
angular resolution is given by the number of angular modes /max and the radial resolution 
by the number of radial modes rimax- We want to keep the resolution constant, so we have 
chosen to link / and nmax(0 such that rimaxlO + V^ = R/rmin] '"min is then the smallest scale 
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linear theory velocity field is easily calculated in terms of these, as 



v(r, c) = Ho(3 ^ 6i^nV f ^mM^^%^) . (2) 

Imn \ '^In ) 

The problem is to determine the ^imn from redshift survey data, especially correcting for 
distortions of redshift space (s, u;) relative to real space (r, w), arising from the velocity field. 
In this section we briefly describe how this is done, in the method introduced by FLHLZ, 
which we extend and apply in this paper. Full derivations are given in the Appendices. 



2.1. Inverting redshift space distortions 

From an all-sky redshift survey (sj, ujj), one can compute sums of the type 

Pfmn = Yl ^is^)Km{^^)jl{klnSi). (3) 

Si<R 

Here w{s) is a weighting function for the galaxies, which we allow to be arbitrary. Now the 
pf^^ seem like Fourier-Bessel expansion coefficients for the density field in redshift space 
(hence superscript S) but they depend on the selection function of the survey — denoted by 
(j){r) — through the sum, as well as on w{s). They are not the expansion coefficients unless 
w{r)(j){r) = 1. In linear theory the pf^^ can be related to the sought-after 6imn by a messy 
but linear relation. 



P n' 



Here 



Oimn = V4:'K / w(t)io{kinr)r dr (5) 



represents a monopole correction: initially we expand the density field. In order to 
transform to the overdensity field, we have to divide by the mean density p and subtract 
the / = (monopole) term from the coefficients. Note that the monopole correction does 
not really depend on m, but we have kept the index for consistency. The second step 
in equation ^ corrects the redshift space overdensity coefficients for the redshift space 
distortions expressed by 



W = J w<PUk^r)UK,rydr 
-(3 / w(t)ji{kinr) 



/(/ + 1) ^\ j[{kn'r) d\n(t)' 



klir"^ I kn'T d\nr 



n 



r'^dr. (6) 



probed by any given mode. This is equivalent to the scheme used by FLHLZ; because the 
zeroes of the Bessel function ji{z) are asymptotically given by zin — 7r(n + //2), their scheme 
of setting a fixed upper limit to z for every ji amounts to keeping nmax(0 + V^ a constant 
for every /. 
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The matrices 

iPi)nn' = J W(t)ji{kinr)ji{kin'r)r'^dr. (7) 

carry information about the weight function w. 

Extracting 5imn from equation @j involves solving a matrix equation for each {1,171). 
Pinn', Oinn', and Zinn' depend on the weight function w{r) and the selection function 0(r) 
but not on data. 

The weight function w{r) could be set to l/(j){r) to eliminate the selection function 
from the formulas. The correction matrices Pinn' would then become diagonal and (0) 
would reduce to the orthogonality relation (|A16|) for spherical Bessel functions. The other 



extreme is to weight all galaxies equally w{r) = 1, leaving the correction of the selection 
effect fully to the matrices Pimn- The difference between these two approaches lies in the 
errors induced by shot noise. The choice w{r) = 1 tends to extrapolate information from 
the well-sampled regions into less-sampled regions, while also propagating shot noise from 
less-sampled regions to well-sampled regions. The choice w{r) = l/0(r) keeps the effect of 
shot noise more local. For the data set we have analyzed in this paper, the difference made 
by w{r) is very small (see figure ^ below) indicating that the errors are nowhere dominated 
by shot noise. However, the situation is different when we encounter an analogous problem 
in Section |^. There the selection function is angular and becomes zero in the unobserved 
region; the information-propagation aspect then becomes crucial for filling in this region. 

The expressions above are equivalent to FLHLZ but not identical because 

• we have chosen to normalise the Fourier Bessel coefficients so that Fimn always have 
the same dimensions as F{r,uj). 

• our coefficients always refer to the Fourier Bessel basis set, whereas in FLHLZ they 
sometimes refer to an intermediate basis set that depends on the galaxy weight w{r). 

• we work entirely in the CMB rather than in the Local Group rest frame. 
For this reason we have given full derivations in Appendix |^. 



2.2. Wiener Filter 

A Fourier-Bessel expansion computed as above will contain spurious extra power from 
shot noise. This spurious power can be suppressed by a Wiener Filter as explained in [Press 
et al (1992)1 , p. 548, where the filter is 



power in signal 
power in (signal + noise) 



The derivation of this expression, however, as given in [Press et al (1992]|, is only vahd 



for a scalar transfer function, whereas in our case the transfer function is Z/nr,'- We will 



therefore have to derive our own Wiener Filter operator. The filter and its derivation are 



given in Appendix M, which mostly follows Zaroubi et al (1995). But we have presented 



the derivation in full since some of our normalizations are different and since we needed to 
preserve the arbitrariness in w{r) throughout. 

The Wiener Filter requires knowledge of the power spectrum P{k) of the underlying 
density fluctuations. We do not know this in general, but since the filter is a correction, a 
first order error in the filter will only introduce a second order error in the full reconstruction 
(cf [Press et al (1992) , p. 548): even a fairly crude approximation of the power spectrum will 



make the filter work. We therefore use a CDM power spectrum with F = 0.4, normalised 
to as = 0.8, as the filtering power spectrum for all our reconstructions. This resembles the 
spectrum fitted to data by Peacock and Dodds (199"4]| and all of the CDM spectra of the 



A^-body simulations (see figure 0). A Wiener Filter on the basis of this power spectrum does 
as well in reconstructing velocities from simulated catalogs as a filter based on the power 
spectrum underlying the simulation in question. The advantage of using this constant filter 
rather than the 'correct' filter based on the underlying power spectrum is that the tests 
on the simulated catalogs will then accurately refiect the error introduced into the PSC2; 
reconstruction since we do not know the correct power spectrum in that case. 

Figure |1| shows power spectra for several cosmologies, including our selected power 
spectrum and the one fitted to existing data by [Peacock and Dodds (1994] . 



3. Mask Correction 

As mentioned above, the formalism developed in the preceding sections is only valid 
for a full-sky catalog. Most redshift surveys are masked in some way, so it is necessary 
to correct for the unobserved parts of the sky. This is done by introducing fake galaxies 
into the unobserved regions while trying to preserve the statistical properties of the galaxy 
distribution of the observed regions. The process is referred to as 'mask correction'. 

The usual method is to clone the fake from the observed galaxies by extrapolating the 
distribution from the observed into the unobserved regions. This works reasonably well, 
but has the following disadvantages: 

• The joining- at-the-seams between the masked and unmasked regions has a tendency 
to introduce spurious power on some scales (since the correlation vanishes at the 
boundaries), which can contaminate the Fourier Bessel coefficients. 

• The cloning method only works if the masked areas are very small compared to the 
unmasked areas. It is hopeless for redshift surveys with very small sky coverage such 
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as the 2dF survey. It is therefore interesting to try out a Fourier Bessel based mask 
correction here for possible apphcation with that survey. Naturally, with small sky 
coverage, only high-/ modes might be usefully constrained, but this may still extract 
information about density and velocity fields on scales smaller than the survey. 



As mentioned in section |A.3| , it is in principle possible to treat the angular window 
function if{u}) at the same time as the selection function 0(r) but it is computationally 
problematic. We therefore choose to correct in two steps: we first treat the mask and then 
the selection function. To this end, we have to first expand the density field in such a way 
as to decouple the angular modes completely from the radial modes. Let g{s,u}) be the 
'raw' redshift space density, i.e., 

g{s,u)=<Pir)p'is), (9) 

separate the angular and radial parts and expand 

0{S, u;) = Y^ gimnYlm{u})jo{konS). (10) 

Imn 

Note that because only jo occurs in this expression, the radial and angular basis functions 
are truly independent and are not eigenfunctions of the Laplacian operator. They do, 
however, form an adequate description of the density field for the purpose of correcting for 
the angular mask. 



We minimise 



V9(l 



E 

Imn 



QlmnYlm{'^)io{kQnS] 



2 



dus'^ds (11) 



to obtain 

I ip{u:)Yi*^jo{konS)du:s'^ds = Cq^ J2 Ql'm'n J <^{i^)Yl'm'{i^)Yi*^{uj)du:. (12) 

ll'mm' 

The integrals can be replaced by sums over galaxies as before. We invert the matrix on the 
right hand side and thereby recover gimn and hence g{s, uj). This recovered density field can 
then be sampled in the masked regions to provide the fake galaxies. 

As an additional way of minimising the errors in this procedure, we invert the matrix 
with a conditioned inversion, i.e., after diagonalising we note all those eigenvalues which 
are less than 1 % of the largest eigenvalue and eliminate the corresponding eigensubspaces. 
The eigenvalues of the matrix will be unity for unmasked modes and small for mostly 
masked modes. Hence, the noise can be suppressed by suppressing the masked modes. The 
procedure is similar to applying a Wiener Filter and ensures that the density field is not 
dominated by spurious features in the masked zones. 

Figure shows a slice of the masked and unmasked galaxies in the x-z plane. 
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Fig. 1. — Power spectra for all models considered: SCDMC (standard CDM with F = 0.5, 
dg = 1.1, n = 1.0), SCDMG (standard CDM, COBE normalised with F = 0.25, ag = 0.55, 
n = 1.0), LCDM (A CDM model with F = 0.25, ag = 0.93, n = 0.3, A = 0.7), and 
PD94 (Peacock & Dodds, 1994), dashed line shows the power spectrum used for our filtering 
(F = 0.4, as = 0.8) 
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Fig. 2. — Projections of galaxy positions onto a plane in x-z (where the x-axis points towards 
/ = 0° and the z-axis towards b = 90°; y e [—20, 20]); unfilled and filled masked zones; Note 
how the mask correction produces voids around (70,0) and (-20,50) which appears to be 
extrapolations of voids in the surrounding data. 
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4. Tests on Simulated Catalogs 

The method described above can now be tested on simulated catalogs with known 
selection function 0(r). We reconstruct the real space density and velocity fields as described 
above and then sample it at the real space galaxy positions. The r^ of each galaxy is 
reconstructed from its Sj and the radial velocity field (cf equation [A23| ) by Newton-Raphson 
iteration. Reconstruction is done for a radius R of 170 /;,~^Mpc and a minimum radial 
resolution r-an-a of 5 /i~^Mpc . Our galaxy weight is given by w{r) = l/(f){r). We will use 
this to be more in line with standard procedure in all of the following reconstructions unless 
otherwise stated. 

Figure |^ shows a comparison of real against predicted velocities in the r, 6, directions 
for an LCDM catalog (for the parameters relating to the model abbreviations LCDM, 
SCDMC, and SCDMG see the caption of p. The solid line is the line of perfect correlation 
and the dashed line represents a fit to the data. Note that this is not a fit in the least 
squares sense. We found that least squares fits have the tendency to be too easily dominated 
by outliers and therefore devised a more stable fitting routine (see also the discussion in 
Press et al (1992]| , p. 700). We constrain the fitted line to pass through the origin and 
choose the slope to be such that if we were to rotate the line by 90° it would cross exactly 
half the points. This does not ignore the outliers but it only accounts for their numbers 
rather than for their distance to the line. 

The figure illustrates that most of the error in the reconstruction lies in the r 
direction, as would be expected since the redshift space distortions reside exclusively in 
that dimension. The error in the 6 and directions is produced by shot noise only. In the 
following we will here only plot the radial velocity comparisons, since they are the most 
interesting in assessing the performance of the method. Note also that the reconstruction 
was performed using a 'constant filter', i.e., instead of using the correct power spectrum for 
this cosmology, we used the chosen power spectrum discussed in section |2.2| . 



Figure ^ shows the effect of using different types of filter. The first panel shows a 





Fig. 3. — Real (vertical) against reconstructed (horizontal) velocities: Vr, vg, v^ for an LCDM 
catalog; plot galaxies out to 60 h~^Mpc , 1/0 weighting, constant filter; it is obvious that 
the maximum error is in Vr- 
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reconstruction for an LCDM catalog using no filter. The second panel shows the same 
catalog, but in this case the velocities were reconstructed with the constant filter. The 
variance is smaller now since most of the shot noise is smoothed away. The third panel 
shows exactly the same case but we have here used a constant weighting function w{r) = 1. 
The result illustrates that in the case of the velocity reconstruction the results are the same 
whether we use an inverse selection function or a constant weighting for each galaxy. 

Figure |^ shows the real against predicted radial velocities for one catalog of each 
cosmology using a constant filter reconstruction. 

Table |^ shows the results of applying our line fitting to the constant filter 
reconstructions. The parameter B is the slope of the line averaged over all 10 catalogs in a 
given cosmology, and d is the average distance of a point from the line. 

A high-density selection in the Fourier Bessel case is more damaging to the correlation 
than a high-radius selection. Note that in the worst possible cosmology (SCDMC), the 
correlation practically ceases in the high density case but is comparatively normal (by the 
standards of that universe) for the high radius case. For all cosmologies, there is a clear 
tendency to overpredict velocities in high density regions. 

Figure ^ again illustrates the dependence of reconstruction error on local density: 
for the same representative reconstructions for each cosmology as in figure |], we plot the 
reconstructed density field (contours) and the difference between the reconstructed and the 
real velocity field (arrows). 

The differences are most marked in the SCDMC case and in all three cases, there is 
a clear correlation of velocity error with overdensity. Since the arrows point towards the 
overdensities, the infall into clusters is overpredicted. 

To test the dependence on radius, we also plot the velocity difference for the SCDMC 
case out to 170 h~^Mpc (figure ^. Note that there is no indication that the greater shot 
noise at high radii produces a systematic effect. Even in those regions, the errors caused 
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Fig. 4. — Real against predicted velocities for three different types of filters; LCDM catalog; 
plot galaxies out to 60 /i^^Mpc , 1/0 weighting in all but the last panel. 



-13- 



-1500 



1500 -1500 



Fig. 5. — Real against predicted velocities for three different cosmologies; one catalog for 
each cosmology; plot galaxies out to 60 h~^Mpc . 



LCDM Results 


Case 




B±aB 


d ± a^ 


r < 120 




0.92 ± 0.08 


244 ± 25 


r < 60 




0.80 ± 0.07 


220 ± 31 


60 < r < 120 




1.14 ± 0.22 


270 ± 35 


6 < 1 




1.17 ± 0.15 


233 ± 25 


6> 1 




0.76 ± 0.34 


288 ± 67 


-1<6 <l,r 


<60 


0.95 ± 0.07 


184 ± 24 


SCDMG Results 


r < 120 




0.90 ± 0.07 


224 ± 26 


r < 60 




0.82 ± 0.07 


212 ± 29 


60 < r < 120 




1.05 ± 0.19 


241 ± 40 


6 < 1 




0.99 ± 0.08 


232 ± 23 


6> 1 




0.74 ± 0.18 


278 ± 52 


-1<6 <l,r 


<60 


0.85 ± 0.07 


193 ± 27 


SCDMC Results 


r < 120 




1.01 ± 0.11 


457 ± 49 


r < 60 




0.84 ± 0.12 


407 ± 71 


60 < r < 120 




1.27 ± 0.13 


517 ± 28 


6 < 1 




1.17 ± 0.16 


467 ± 28 


6> I 




1.37 ± 1.51 


751 ± 685 


-1 < 6 <1, r 


<60 


0.92 ± 0.15 


372 ± 77 



Table 1: Parameters of fitted line and velocity dispersions for three different cosmologies 
and 6 different cases. All distances are in h'^Mpc . B is the slope of the fitted line, as the 
dispersion in the average over all 10 catalogs, d is the average distance of a point from the 
fitted line and ad its dispersion. 
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Fig. 6. — Av = fpred — "Wreai for three different cosmofogies (same catafogs as before); Density 
levels indicate overdense regions only, bold line is at 5 = and the contours are evenly 
spaced in 13 steps from S = 1 to 6 = 7. Velocity vectors are drawn on an arbitrary scale. 
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by overdense regions far outweigh any radius dependence. Again, this is a trend aheady 
indicated in table Bl. 



5. Local Group Velocity and (3 

The velocity field reconstructed by the above method is of particular interest at the 
origin: by comparing the reconstructed velocity, v(0), of the Local Group (LG) to the Local 
Group velocity, wlgi measured from the dipole in the CMB we can recover (3 . 

Since only the dipole contributes to the velocity v(0) (cf equation |A3CI| ), we can set 
^max = 1 for this particular calculation. We also increase rmin to 10 /i~^Mpc to increase 
the stability of the reconstructed LG velocity. If the resolution is too high, the calculated 
velocity tends to be too easily infiuenced by small but nearby density fiuctuations. 

Figure || shows the measured LG velocity amplitudes for reconstructions using /? = 
0.3 to 1.0 (steps of 0.1) rescaled to /? = 1. We plot this against the angle between the 
reconstructed LG velocity direction and the CMB dipole direction (misalignment angle). 

Note that the reconstructed velocity does not seem overly sensitive to the (3 assumed 
in the reconstruction. The rescaled amplitudes are all in the range between 1000 kms~^and 
1100 kms~^and hence we can conclude that (3 ~ 0.7. The very good alignment between the 
reconstructed velocity and the CMB dipole indicates that we are satisfactorily sampling the 
matter that causes the acceleration. Note also that the convergence of the dipole amplitude 
is very good. The Fourier Bessel method works entirely in the CMB frame and we therefore 
do not have any signature from the Kaiser rocket effect (discussed, for example, in |Straus£| 
et al (1992)1 ). 



Instead of performing a likelihood analysis, in this case we evaluate systematic effects 
by calculating the LG velocity for all 30 catalogs to look at the scatter in 

(v(O)-vlg)'. (13) 

We render this dimensionless and define 

/^'=(f— ) -2cos?9— + 1), (14) 

\\vlgJ vlg ) 

where °d is the misalignment angle. This method may be less sophisticated than a likelihood 



analysis (e.g., [Schmoldt, et al (1998)|) but it makes fewer assumptions about the underlying 



cosmology. It works on the simple basis of comparing real data against results from 
simulated catalogs. 

Figure || plots the quantity v/vlg against the misalignment angle for all simulated 
catalogs and the PSC^; together with contour lines of the surface equation 111. The LG 
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Fig. 7. — Av = f pred — "^reai for SCDMC; density levels indicate overdense regions only, bold 
line is at (5 = 0, thin lines are overdensities from 5 = 2 to 5 = 32 in 16 steps; the smallest 
velocity difference for this plane is 18 kms~^and the largest 1788 knis^^. 
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Fig. 8. — Dipole amplitude and misalignment against different % 
the velocity is rescaled to /9 = 1 in the upper panel. 



for the PSC2;; note that 
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velocity for the PSC2 was reconstructed using (5 = 0.7 as was indicated by the above result. 
The inner contour encloses 2/3 of the simulation points and hence it defines our 68 % 
confidence limits on f3 . Therefore, we recover f3 = 0.7 ± 0.5. This result seems a lot less 



restrictive than the one quoted in Schmoldt, et al (1998) but it has to be noted that in that 



work the results separate the high and low normalisation cases. In this case, we consider all 
simulations at the same time - drawing together the two different values for j3 in ^chmoldt 
et al (1998)1 would give a similarly ill-constrained result. 



It is interesting to note that most of the velocity ratios are close to unity but that 
the misalignment angle can get quite large. This is particularly the case for the SCDMC 
catalogs, which are (as already mentioned) dominated by large fluctuations at small radii. 
This seems to affect the misalignment much more than it does the amplitude. 



6. Velocity Anisotropy and P 

The density field 6^^^ should in general be perfectly isotropic, i.e., there should be 
nothing special about the radial direction. The radial redshift space distortions in the 
density field 6^ , however, depend on (3 so that j3 can in principle be recovered from analysing 
those distortions. 

We recover the full 5^^ in our method, and if we do so with the correct f3 , the density 
field should be isotropic. Therefore, we simply devise an anisotropy test which will detect 
radial anisotropics caused by a reconstruction of the density field using the wrong /3 . Note 
that this method tests reconstructed real space velocities and not the redshift space density 
field. In order to avoid problems caused by a small sample (i.e., in order to improve the 
chance that the normal anisotropics will average out in the volume considered) we restrict 
the analysis to regions where 6 < 1. 

Increasing /? in Zinn' (in equation |^) suppresses structure in the radial direction (since 
the redshift space distortions create such structure artificially). If our assumed (3 is too 
high, then the recovered Vr will be too low. We therefore define the parameter 

3 {V^} - {Vr}^ 

where the averages are weighted by the local density (1 + 6). We expect r] to pass through 
unity at the correct /3 . 

To calculate the value of 77 for different catalog reconstructions, we sample the velocity 
field in a radius of 60 h~^Mpc of the origin with a collection of 800 random points having 
6 < 1. (We exclude regions with 6 > 1 since our reconstructions have larger errors there.) In 
figure |l^, we plot r] against various trial /3 for 10 SCDMG and 10 LCDM catalogs. We have 
discounted the SCDMC catalogs for this test because of their unrealistically large density 
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Fig. 9. — Dipole determinations for all catalogs: triangle = PSCz, + = SCDMC, 
SCDMG, square = LCDM. 
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fluctuations on small scales (much larger than observed clustering) and partly because of 
the high variance in the predicted velocities (see table 1). An anisotropy test on SCDMC 
would be too noisy to give useful results. 

In general, the scatter is too large for a clear statement about the j3 at which the lines 
cross unity. However, we can recognise some trends. Generally, the SCDMG points (/? = 1) 
are lower, i.e., they cross unity later than the LCDM points {(3 = 0.5) as would be expected. 
This suggests treating 77 as a statistic, and the values obtained from the 10 SCDMG and 10 
LCDM catalogs as samples of the distribution of rj in the respective cosmologies. We find 
(cf p!0|) that, for all trial (3 , the value of 77 from PSC2; exceeds at least 9 (usually all 10) of 
the values from PSC2;. This is not the case for LCDM. We therefore conclude that SCDM 
is excluded at the 90% confidence level. 

It is worth noting that not all (3 =1 models can be rejected on the basis of this result, 
since models with more large-scale power will also have more cosmic variance, increasing 
the scatter in rj. 



7. Local Density and Peculiar Velocity Fields 



Figures |ll| and show the PSCz velocity and density field in the supergalactic plane 
out to 70 and 170 h~^Mpc respectively. The general features agree very well with the 
velocity field reconstructed in Pranchini et al (1998)| . We recover the same continuous flow 



from the northern end of the Perseus-Pisces cluster to the GA region. However, we also find 



a back-infall into this region, which was not observed in Branchini et al (1998) . This is more 
apparent in figure |T^: there is a clear division between the back-infall into Hydra-Centaurus 
and the subsequent infall into the Shapley cluster (-100,90). The highest velocity for any 
of the back-infalling galaxies is 494 kms~^and therefore above the dispersion level given 
in table |l|. However, table |l| also indicates that the infall into high density regions will 
in general be overpredicted. Likewise, comparing the real and predicted velocity maps 
of the simulated catalogs, we find that the overpredicted infall makes for a more sharply 
peaked velocity field in the reconstructions and hence in about half of the maps for the 
SCDMG catalogs, for example, it is easy to identify regions where the reconstructions show 
a back-infall into a cluster that does not exist in the real velocity field. We cannot therefore 
consider the evidence conclusive. 



Figure |T3| presents the PSCz density field. We have plotted all those galaxies, for 
which the local density S{r) as determined from the density coefficients is higher than the 
threshold 6 = 1.0 and which are within a radius of 130 h~^Mpc of the observer. 

We labelled all the known superclusters as before. The same structures are identified 
and it is particularly interesting to note the extension to Hydra. 



p 
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Fig. 10.— Plot of T] for different trial (3 from PSCz, and from 10 SCDM and 10 LCDM 
catalogs. The PSCz line lies among the LCDM values but is almost completely disjoint from 
the SCDMG values. 
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Fig. 11. — Velocity field of the PSC^ survey: radius out to 70 h^^Mpc . Contour lines 
denote density field (fat line = zero density contrast, broken lines = underdensities, thin 
solid lines = overdensities) . 
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Fig. 12. — Velocity field of the PSCz survey: radius out to 170 h ^Mpc 
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Fig. 13. — Clusters in the PSC2; survey: Density level at 1.0. 
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8. Conclusions 

We have corrected the redshift space distortions in the PSC2; survey via the Fourier- 
Bessel method of FLHLZ with two extensions: a generic galaxy weight and a matrix 
correction for the masked zones. The method was tested extensively on mock catalogs 
extracted from A^-body simulations. It was shown not to have any systematic biases except 
for an overprediction of velocities in high density regions. We find 

• an LG velocity that well reproduces the CMB dipole for /? = 0.7. We estimate the 
error on this value from the scatter in the simulations as ± 0.5. 

• that an Q = 1., cluster normalised and standard CDM cosmology is ruled out with a 
confidence of 90 % from a new anisotropy test. 



a flow field in the supergalactic plane that is consistent with the results of [Branchim 



et al (1998)1 apart from an observed back-infall into the Hydra- Cent aurus region. 
However, we treat this result with caution since the simulations show spurious features 
of this type. 

a density field that shows all the usual clusters and confirms two of the ones identified 
by |Webster, Lahav, and .b'isher (1997]|. There is a possible conincidence of Saar2, 



a structure identified in the PSC2; by [Saar (1996)| , and the C( cluster of [Webster, 



Lahav, and Fisher (19971 and likewise of Saarl with the Orion and A539 clusters. We 
also confirm the existence of A3627 and identify an extension to the Hydra-Centaurus 
supercluster. 



I. M.S. and P.S. wish to dedicate this paper to the memory of Veikko Saar who started 
this work but was unable to finish it. The authors are grateful to Bill Ballinger, James 
Binney, and Radek Stompor for helpful feedback. I. M.S. acknowledges financial support 
from Oriel College and PPARC, and is particularly grateful to the Sasakawa Foundation 
and Tokyo University for funds and hospitality. 



A. Dynamical Theory 

This appendix derives the equations in section ^ that express the express the real-space 
overdensity and velocity fields under linear theory from survey data via a Fourier Bessel 
expansion. 

We will first introduce the mathematics associated with the distortion correction. 
These are independent of the choice of basis functions so it is only after the derivation 
that we motivate our particular set of functions (spherical harmonics and spherical Bessel 
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functions) by the observation that they are eigenf unctions of the Laplacian operator and 
therefore render the velocity calculations particularly easy. We then present the form of the 
distortion correction for this basis function set. 



A.l. Redshift Space Distortions 

The transformation between redshift space s and real space r is given by 

s = r + u;u{r), (Al) 

where u{r) is the radial peculiar velocity at a point r and w is the angular unit vector, i.e., 
uj = r. 

Consider some function /(s) of the galaxy positions s in redshift space given by a 
survey. Summing over that function at the galaxy positions will be equal to the integral 
over all space of the function multiplied with the the sampling function of the galaxies, i.e., 

Y: /(s.) = / rfa; rf{s,u.)p'{s,u.)<f>{r)s'ds, (A2) 

where p^{s,uj) is the redshift space density field, 0(r) is the selection function (note that 
it depends on r, not s, since a galaxy's inclusion into a fiux-limited survey depends on its 
position in real space) and R is the maximum radius of the survey in redshift space. To 
a zeroth order approximation, this sum will be equal in redshift and in real space. In the 
following, we will work out a first order correction to that approximation by finding the 
analogous sum in real space. 

We now hold a; constant since the effect of redshift space distortions is isolated in the 
radial component. Mass conservation requires that 

p^{s)s^ds = p^'irydr, (A3) 

where p^*^ is the real space density field, so we can rewrite the right hand side of equation 

f{s)(f){r)p'^{s)s^ds= f "^ f{s)(j){r)p^^{rydr. (A4) 

Jo 

Note that the boundaries on the integrals have changed. The maximum radius R is defined 

in redshift space and therefore also has to be transformed to real space. We then expand 

f{s) in a Taylor series of the form 



fir + u{r)) = f{r)+uf'ir) + ... (A5) 

and drop all second order terms, so that 

R pR—u pR—u 

f{s)p^{s)cf>{r)s^ds = / f{r)p^^{r)^{rydr + / uf'{r)p^^{r)^{rydr. (A6) 
Jo Jo 
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The first integral on the right hand side has to be spht into two parts 



R~u 



/(Pp'^^r'dr 



Jr^ 



(A7) 



since it is a zeroth-order term and the integral between R and {R — u) cannot be neglected 
(note that we have ceased to give the dependence on r in all the functions to make the 
equations neater). The last term in equation [A7| is first-order, so we can replace p^'^ by its 
average value p, multiply the integrand with u, and evaluate the functions at R. 

The second integral in equation ^^ is a first-order term, so we can again replace p^*^ 
by p and integrate by parts. Simplifying, we get 



f{s)p''{s)ct){r)s''ds 



R f^ d / \ 

/(pp^'^r'^dr - p f—- (u(f)r^) dr 
Jo dr ^ ' 



(A8) 



Here the first term is the real space expression (analogue of the redshift space sum in 
equation |A2[ ) that we were looking for and the second term represents the redshift space 
distortions. We can gain some insight into the nature of these distortions by rewriting this 

term as 

' d\a.0 \u 

M' + ■ 



P 



R 



d\a.r 



2 - 

/ r 



dr. 



(A9) 



where the terms in the square brackets now describe the various contributions to the 



distortions, i.e.. 



the galaxies' peculiar velocities: the velocity field at a real space position r will be 
different from the field at s - hence the m' term. 

the selection function: there is a change in the selection function between r and s - 
hence the (iln0/c?lnr term. 

the volume change: a shell at r has a different volume from a shell at s - hence the 
2ulr term. 



These were first pointed out by [Kaiser (1987)| . Note that of the two kinds of redshift space 
distortion discussed in that paper, this expression only treats the second, linear kind and 
hence a redshift space correction on the basis of this expression will only work in the linear 
regime. 



Now, putting back the a;- dependence, we have 

^ /(s,) = \d^ r/(r,u;)0(r)p«^(r,u,)r2rfr- 

P fir,u:)— {u{r, u)(j){r)r^j dr. 



Si<R 



(MO) 
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This represents the general equation for the correction of the redshift space distortions. 
We shall find a more specific expression after deciding on the form of /. 



A. 2. Peculiar Velocities 

In linear theory, the velocity is the gradient of the gravitational potential (see Peebles| 



T993)i p. 116) 

v = -^f^V^. (All) 

Note that v is, as usual, a comoving peculiar velocity at the present time of which u in 



equation ^ is the radial part. Poisson's equation states that 



VV = ^no5ra.tterHl (Al2) 

We use the standard approximation f{Qo) ~ ^o^, convert the density contrast (^matter to the 
galaxy density contrast 6, assuming constant bias. The velocity therefore becomes 

V = Ho(3V{V-^6), (A13) 



A. 3. Fourier Bessel Expansion 

Since peculiar velocities are our main interest, we want to render the solution of 
equation |A13| as easy as possible. It is obvious that the calculation of v will be trivial if 
we can find an expansion for 6 that is an eigenfunction of the Laplacian operator. We 
also require the expansion to be in spherical coordinates]^ : firstly because a separation 
of the angular and radial parts will concentrate the redshift space distortions in only one 
dimension, but also because the decreasing selection function of a redshift survey will tend 
to render its volume spherical. We therefore choose the Fourier Bessel expansion 

F(r,a;) =J2FUk)YiU^)Jji{kr)dk. (AM) 

Im 

Since the volume of the redshift survey is finite, we approximate the integral on the right 
hand side by the sum 

F(r, uj) = Yl FlmnYlrniuj)jlikinr), (Al5) 

Imn 



^^Spherical coordinates are given by x = r sin 6' cos = r cos 6 cos/, y = rsin^sin0 = 
r cos b sin I, z = r cos 6 = r sin b. We continue to use a; as a shorthand for 6,(j) - (p here being 
the polar angle, not the selection function. 
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i.e., we confine ourselves to A;- values satisfying certain boundary conditions on the survey 
range R. There will then be an orthogonality relation 

r 3i{hnr)ji{hn'r)rHr = 5f„,Q;\ (A16) 

where 5^ denotes a Kronecker 5. The Yim{ijj) are of course orthonormal. The values of 
Cin and kin depend on the precise boundary conditions. We follow one common choice, 
which is to require 5 = at r > /? with 5 allowed to be discontinous at r = i?, but the 
logarithmic derivative of the gravitational potential dhiil) / d\nr required to be continous at 
that boundary. In this case, 

kin are zeroes oi ji^iikR) (A17) 

Crn' = \R\3l{hnR)f. (A18) 

For a discussion of other possible boundary conditions see FLHLZ, appendix A. 

If F{t) was known everywhere, we could simply invert equation |A15| using orthogonality 
to obtain 

Fimn = Cm I Y;„^nMi{kinr)F{r)dh (A19) 

JVh 

However, F{r, 6, (p) is known only with some position-dependent accuracy W{r) (such as, 
for example, the selection function combined with an angular mask). We then estimate 
Fimn in a least-squares sense by minimising 

2 

-^(r) - Yl ^l'm'n'Yl>„i>{u})ji>{ki>n'r) 
^« I'm'n' 

Differentiating with respect to the -F^^n and setting the result to zero, we obtain 

/ F{r)Yi*„ji{kinr)W{r)dh= ^ F.^'n' f W{r)Y,^>Y;^iik,n'r)3i{kinr)dh. (A21) 

This is the most general possible expression. We can recover the Fimn (corrected for the 
function W{r)) by calculating and inverting the matrix on the right hand side. However, 
since the angular part will have (2/max + 1)^ elements and the radial part nmax(O) the matrix 
will become huge and computationally difficult to invert. This problem can be alleviated 
if W depends on r only. In the following, W{r) = (f){r)w{r), where w{r) is some weight 
function. In that case, the integral over YimYif^, will reduce to Kronecker deltas because 
of the orthogonality condition (spherical harmonics are orthonormal on a sphere) and we 
obtain 

F{r)Y*ji{kinr)(p{r)w{r)dh = Y,Fin,n' / 0(r)«;(r)j;(A;;„,r)ji(A;;„r)r2rfr. (A22) 

We can recover the Fi^n by multiplying by P/nn' on the left hand side, where the inverse of 
Finn' is defined in (^. This will enable us to correct for the selection function but not for 
any angular mask. Mask treatment is discussed in section 0. 



W{r)dh. (A20) 
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A.4. Application of Expansion 



We now use the Fourier Bessel expansion (|l]) and calculate the version of equation |A10 



valid for this choice of basis functions. Note that equation |l| is just a version of equation 
A15| with the density contrast field 6 as the function F. Equation |l] together with equation 



A13| leads to 

uir, u) = Ho/3 Y: ^-^^j^YU-^)5i^n, (A23) 

Imn ™ 

and 

^ + — = ifo/5^ f^T^ - ll Jliknr)Y,m{-^)S,mn. (A24) 

dr r ;^V kLr^ J 

For the last step, we have used the Bessel equation to eliminate i"{kinr). 

We now choose the function /(r, w) in equation [A1U| as 

f{r,u:)=w{r)ji{h^r)Ya-^), (A25) 

where w{r) is a weight function (e.g. 1 or l/0(r) or anything else as discussed above). 
Substituting p = p{l + 5) with 5 given by equation || into equation |A10| , we obtain 



- ^ w{si)ii{kinSi)Yim{u}i) - V47r / wcpjoikinr-y'^dr 
Y.^i"^n' / W(t>ji{hnr)ji{kin'r)r'^dr 



-/?E^w/-0jKA:.nr) 






r^dr. (A26) 



Here we have moved the / = terms to the left hand side. In a shorter notation, this 
becomes equation (^. 



A. 5. The Local Group Velocity 

The procedure outlined above will allow the reconstruction of the entire velocity field, 
but of particular interest to us is the velocity of the Local Group. However, our methodology 
uses polar coordinates which are not defined at the origin. To solve this problem, we could 
either simply use points close to the origin or evaluate the hmit by explicitly adding up the 
contributions to the acceleration from the surrounding density field. We have chosen the 
latter method and our treatment follows closely that of FLHLZ. 

The LG peculiar velocity is given by 
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Substituting a Fourier Bessel expansion for the density field in the usual way (cf equation 
A15D , we obtain 



v(0) = T^ E t n^ dr'jiihy) j dJJYiU^'), (A28) 

47r r-^ JO J47r 



Imn "" ■^^'^ 



where rmax is the radius out to which the density field is considered for the calculation of 
the LG velocity. We can evaluate the last integral by noting that 

u,- f dJJYim{u^') = ^6^, Y, F;™M, (A29) 

"'4'r O m=-l,0,l 

(see FLHLZ) where S^ denotes the Kronecker delta. Hence, only the dipole term / = 1 
survives. Our choice of boundary conditions implies kin = ^, and we can do the first 
integral analytically to obtain 

v(0) = ^=E(-v^3fJ(5n„)i+v^S(5n„)y + 3fJ(,5io„)z) 

Afi-^sinf!^^^^)) (A30) 

where 3? and Q" refer to the real and imaginary parts of a complex number respectively. 



A. 6. Computational Comments 

In the code, we calculate the Yi^ and jin in the standard way using Numerical Recipes 
( Press et al (1992)1) routines. The first derivative of the Yim can be calculated using 
recursion relations from [Abramowitz and Stegun (1965)| , equation 8.5.4. Since 11, -m = Y^^ 



we need to store the coefficients Fi^n only for m > 0. 



B. The Wiener Filter 

As shown in Appendix 0, the real coefficients 5^^ ^^^ the measured coefficients, say 
Cimri) are in principle connected in a relationship of the type 

Clmn = / .Zinn'SjJ;„,, (Bl) 

n' 

where Cimn is a shorthand for the left hand side of equation ^. However, in the real case, we 
still have the shot noise to deal with, i.e., the relationship above has to be amended to 

Clmn = 2J Zlnn' Ulmn' + f^'lmn') , (B2) 
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where fiimn are the Fourier Bessel coefficients of tlie sliot noise. 

As a first step, we can simpfify tilings by noting that there is no coupling between 
different {l,m). Therefore, for a given {l,m), we can write equation IB^ more concisely as 



C=Z(^ + m), (B3) 

where C, S, M are column vectors of size nmax(0 and Z is a nmax(0 ^ "^maxlO matrix (we have 
dropped the indices / for clarity). 

We now want to derive an operator T which will allow us to get a good estimate S^st of 
the real S 

^est = TC. (B4) 

As always, we minimise the difference between the ^est and the real S, i.e., we minimise 

((^-TC)^(^-TC)). (B5) 

In this and the rest of this section, angular brackets indicate an average over different 
realisations of the noise instead of the usual spatial average. Differentiating with respect to 
T and setting the result to zero, we obtain 

T=(^C^)(CC^)-^ (B6) 

We know that (/x) = (/x"'') = since the noise averages to zero, and so, from equation 



B3| we have 

{5C^) = {5S^)Z^ (B7) 

and 

(CCn = Z [(Mt) + (;.;.t)] Zt. (B8) 

We call {SS'^) the signal matrix S (recognising the normal correlation function) and (m/x^) 
the noise matrix N. Hence 



(B9) 



The signal and noise matrices are given by 

S + N=((^ + /x)t(5 + ;x)) (BIO) 

We return to index notation and therefore obtain 

Slmnn' + Nimnn' = i^lmn^lmn') ^ (BH) 





T 


= SZ^ [Z (S + N) Z^" 
= S(S + N)-'Z-\ 


which 


is the form of the Wiener Filter in (H) . 
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where Sj^^^ includes shot noise. Now, consider our estimate of the real space density 
coefficients 



Hmn 



Vr 



<P{r)w{r)5''\r,^)3i{hn)YiU^)Pin. 



(B12) 



We also know that the expectation value for density fluctuations with the shot noise taken 
into account is given by (cf [Bertschinger (1992)| ) 



iRe 



{5^\r,)5^\v-,))=i{\vr-v-,\) 



p(j){r] 



^g^ri-r^) 



(B13) 



^(3) 



where S^ describes a three-dimensional Dirac delta function. The first term in this 
expression describes the signal, whereas the second part describes the noise. ^(| ri — r2 |) is 
the correlation function of the density field. It depends on distance only, since we assume 



isotropy of direction. Note that the cross-term in equation pi3| was dropped since we 
assume that signal and noise are uncorrelated. 



Using equation |B12| to calculate the expectation value equation |B11| and substituting 



equation B13, we can isolate the signal and noise matrices as 



^(1 Tl — r2 \)PlninPln2n' 



Nlmnn' = Y. drr^4>W^Jl{knir)jl{kin2r)PlninPln2n'- 

n-ino •'^R 



(B14) 
(B15) 



The noise matrix can then be worked out explicitly. It does not depend on m. 

The signal matrix can be simplified further by first noting that it cannot depend on 
(f) or w{r) since it describes the signal only. We are therefore free to choose any weight 
function without loss of generality and set w{r)(f){r) = 1 in the following. Additionally, 
consider the relation between the correlation function ^(| ri — r2 |) and the power spectrum 
Pik) 

e(l ri - r2 I) = -^ J rf3kP(A:)e-'^-(^^-^^) (B16) 

and the Rayleigh expansion of the exponential in spherical waves (cf [Arfken (1985)| , p 665) 

e^k-r ^ 47rJ2^ij^(kr)Y,l{u)YU-^,). (B17) 



This yields 



S, 



Imnn' 



Im 



2 r°° r^ 

- / dkk^P{k) / drirlji{kinri)ji{kri) 

TX Jo Jo 

R 

dr2rlji{kin'r2)ji{kr2)CinCin' 





(B18) 
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since the spherical harmonics are orthogonal. The last integral on the right hand side can 
be approximated when the upper limit is set to oo (i.e., we assume that the volume inside 
R is representative of the universe as a whole), and hence 



dr2rlji{kin'r2)ii{kr2)CinCin' 



TV 1 



2k 



^5\){k-kini) 



We therefore have 

which can then be reduced to 
and again does not depend on m. 



~ - / dkk^P{k) / drirlji{kinri)ji{kri] 
n Jo Jo 



n 1 
2fc2 



sh{k-k 



In' 



s, 



Imnn' 



^lnJ^^\kn 



(B19) 

(B20) 
(B21) 
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